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Abstract 


This article describes a search for high-mass resonances decaying to a pair of photons using 
a sample of 20.3 fb”* of pp collisions at = 8 TeV recorded with the ATLAS detector at 
the Large Hadron Collider. The data are found to be in agreement with the Standard Model 
prediction, and limits are reported in the framework of the Randall-S undrum model. This 
theory leads to the prediction of graviton states, the lightest of which could be observed at 
the Large Hadron Collider. A lower limit of 2.66 (1.41) TeV at 95% confidence level is set 
on the mass of the lightest graviton for couplings of k/Mpi = 0.1 (0.01). 
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1 Introduction 


One of the primary goals of the experiments at the Large Hadron Collider (LHC) is the search for new 
phenomena which might manifest themselves in the high-energy regime made accessible for the first 
time by this machine. This article reports on a search for new high-mass resonances in the diphoton 
channel using the full data set of proton-proton (pp) collisions recorded with the ATLAS detector [1] at 
- 8 TeV in 2012. 

The diphoton channel is important for searches for physics beyond the Standard Model (SM) since new 
high-mass states decaying to two photons are predicted in many extensions of the SM, and since searches 
in the diphoton channel benefit from a clean experimental signature: excellent mass resolution and mod¬ 
est backgrounds. One popular example of theories predicting a new high-mass diphoton resonance is 
the Randall-Sundrum (RS) model [2]. This model attempts to solve the so-called hierarchy problem of 
the SM, i.e. the fact that the electroweak scale is 16 orders of magnitude smaller than the Planck scale 
(Mpi w 10^^ TeV). The hierarchy problem poses the question of the naturalness of the SM, and despite 
significant efforts over the past few decades, this problem has not been solved so far. In the RS model, the 
observed hierarchy of scales is accommodated assuming an extra spatial dimension. In this paradigm, the 
fundamental strength of gravity is comparable to the strength of the electroweak interaction, but gravity 
appears much weaker than the other interactions because it is diluted by the presence of the extra spatial 
dimension. 

The RS model postulates a “warped” five-dimensional spacetime configuration. In this configuration there 
are two 3-branes that are the boundaries of the five-dimensional spacetime. The 3-branes can support 
(3 -I- l)-dimensional field theories. The SM fields are locafed on one of the two 3-branes, the so-called 
TeV brane, while gravity originates from the other brane, called the Planck brane. Gravitons propagate in 
the five-dimensional bulk. The fifth dimension is compactified with length r^. The solution to Einstein’s 
equations for this configuration is given by a spacetime metric that couples the standard four dimensions to 
the extra dimension with an exponential “warp factor”. In this configuration, a field with the fundamental 
mass parameter tuq will appear to have the physical mass m = mo, where k is the curvature scale 
of the extra dimension. Scales of the order of 1 TeV are therefore generated from fundamental scales 
of the order of Mpi if kr^ ~ 12 [3]. The compactification of the extra dimension gives rise to a Kaluza- 
Klein (KK) [4] tower of graviton excitations G* - a set of particles on the TeV brane with increasing mass. 
The phenomenology can be described in terms of the mass of the lightest KK graviton excitation {mo*) and 
the dimensionless coupling to the SM fields, k/Mpi, where Mpi = Mpi/ is fhe reduced Planck scale. 
The masses and couplings of the individual graviton states are determined by the scale A;^ = ~ 

(9(TeV) [3]. Values of kjMpi in the range [0.01,0.1] are preferred from theoretical arguments [3]. The 
lightest graviton state is expected to be a fairly narrow resonance for kfMpi < 0.3 [3]. Its natural width 
varies like the square of kjMpi. For klMpi = 0.1, the natural width increases from ~ 8 GeV at mo* = 
800 GeV to ~ 30 GeV at mo* - 2200 GeV. In the first example, the natural width and the experimental 
mass resolution are similar, in the second example the natural width dominates the total width ^ of the 
mass peak. 

Searches for high-mass diphoton resonances at the LHC have been reported by both the ATLAS and CMS 
collaborations, and the latest results are summarized in Table 1. The most stringent experimental limits 
on RS gravitons to date have been obtained at the LHC, using a combination of the dielectron and dimuon 
channels. Limits from searches for narrow resonances in other channels are available in, for example, the 


* The experimental mass resolution at mc> = 2200 GeV is only marginally (6%) better than at mo* = 800 GeV. 
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Table 1: Overview of previous limits on the mass of the lightest RS graviton state, expressed at 95% confidence 
level. The notation {*(~ refers to a combination of the dielectron and dimuon channels. The CDF and DO experi¬ 
ments operated at the Tevatron pp collider. 


Experiment 

[TeV] 

Luminosity [fb 

Channel 

Limit 

Reference 

ATLAS 

7 

4.9 

yy 

2.06 (1.00) TeV for yf/Mpi = 0.1 (0.01) 

[5] 

CMS 

7 

2.2 

yy 

1.84 (0.86) TeV for kjMn = 0.1 (0.01) 

[6] 

ATLAS 

8 

20.5 

ur 

2.68 (1.25) TeV for kjMn = 0.1 (0.01) 

[7] 

CMS 

7 and 8 

5.3 and 4.1 

ur 

2.39 TeV for fc/Mpi = 0.1 

[8] 

CMS 

8 

4.0 

Jet-jet 

1.45 TeV for ;t/Mpi = 0.1 

[9] 

ATLAS 

7 

4.7 

VL+iy- 

1.23 TeV for ;t/Mpi = 0.1 

[10] 

CDF 

1.96 

5.7 and 5.4 

and yy 

1.06 TeV for fc/Mpi = 0.1 

[11] 

DO 

1.96 

5.4 

e^e~ and yy 

1.05 (0.56) TeV for yf/Mp, = 0.1 (0.01) 

[12] 


dijet and WW channels, but these channels are not competitive for the specific case of RS gravitons. The 
limits obtained using these channels are included in Table 1, along with those from earlier searches at the 
Tevatron. 


2 The ATLAS detector 


The ATLAS detector is described in detail in Ref. [1]. It consists of an inner detector (ID) surrounded 
by a solenoid that produces a 2 T magnetic field, elecfromagnefic (EM) and hadronic calorimefers, and a 
muon specfromefer which employs air-core foroidal magnefs. In fhe pseudorapidify ^ region \t]\ < 2.5, fhe 
ID provides measuremenfs of fhe fracks originating from fhe primary pp collision as well as fracks from 
secondary verfices, permiffing fhe efficienf identification of phofon conversions. The EM calorimefer 
(ECAE) is a lead/liquid-argon (EAr) sampling calorimefer. If is divided info a barrel section covering 
fhe region I 77 I < 1.475 and fwo endcap sections covering 1.375 < \t]\ < 3.2. The ECAE is segmenfed 
longifudinally in shower depfh info fhree layers for \r]\ < 2.5 and fwo for 2.5 < \r]\ < 3.2. Up fo \r]\ < 2.4, 
fhe firsl layer consisfs of highly granular sfrips segmenfed in fhe rj direction for efficienf candidafe-by- 
candidafe discriminafion befween single phofon showers and fwo overlapping showers originafing from 
TT*’ decay. The second layer, wifh a fypical granularify of Ap x A(p = 0.025 x 0.025, collecfs mosf of 
fhe energy from phofon showers. Significanf energy deposifed in fhe fhird layer is an indicafion for 
leakage beyond fhe ECAE from a high-energy shower, and fhe measuremenfs from fhe fhird layer are 
used fo correcf for fhis effecf. A fhin presampler layer, insfalled in fronf of fhe ECAE and covering 
fhe pseudorapidify inferval \p\ < 1.8, is used fo correcf for energy loss before fhe ECAE. In fhe region 
\p\ < 1.7, fhe hadronic sampling calorimefer uses iron absorbers and plasfic scinfillafor files as fhe acfive 
material. In fhe region 1.5 < I 77 I < 4.9, EAr is used as fhe acfive maferial, wifh copper or/and fungsfen 
absorbers. The muon specfromefer feafures fhree sfafions of precision chambers fhaf allow for accurate 
measuremenfs of fhe muon frack curvafure in fhe region I 77 I < 2.7. Easf defectors enable muon friggering in 
fhe range \p\ < 2.4. A fhree-level frigger sysfem is used fo selecf fhe evenfs fo be recorded for subsequenf 
offline analysis. 


^ ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point in the center of the detector 
and the z-axis along the beam pipe. The x-axis points from the interaction point to the center of the LHC ring, and the ^-axis 
points upward. Cylindrical coordinates (r, (p) are used in the transverse plane, (p being the azimuthal angle around the beam 
pipe. The pseudorapidity is defined in terms of the polar angle 6 as ;/ = - In tan(0/2). 
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3 Event reconstruction and selection 


The events used in this analysis were recorded using a diphoton trigger. At the highest, software-based 
level of the trigger system, this trigger requires at least one photon with transverse energy > 35 GeV 
and at least one other photon with > 25 GeV. Both photons need to satisfy requirements on the shape 
of the energy deposit in the calorimeter. This trigger is nearly fully efficient for diphoton pairs that pass 
the offline event selection. 

In the offline analysis, only events that pass the standard ATLAS data quality requirements are considered. 
Photons are reconstructed from fixed-size clusters of cells in the EM calorimeter. To ensure good dipho¬ 
ton mass resolution, photons are required to be located inside the precision region of the EM calorimeter 
(|? 7 l < 1.37 or 1.52 < \t]\ < 2.37). They are further required to satisfy > 50 GeV as well as the loose 
identification requirements from Ref. [13] updated for 2012 run conditions. The two highest-E!^^ photon 
candidates that pass these requirements are retained for further analysis. Events are retained for the main 
analysis if both candidates also satisfy the tight identification requirements [13] updated for the 2012 run 
conditions. These requirements include upper limits on the energy leakage into the hadronic calorimeter 
and on the width of the shower in the first and second layers of the EM calorimeter. Events where at 
least one candidate fails the tight photon identification are retained for the study of backgrounds with 
misidentifled (fake) photons in data control samples. 

The invariant mass of the diphoton system is evaluated using the photon energies measured in the calorime¬ 
ter as well as the direction of the photon momenta determined from the positions of the photons in the 
calorimeter and the position of the primary vertex of the hard interaction. Primary vertices are recon¬ 
structed from the tracks in the ID and are required to be formed from at least three tracks. The primary 
vertex of the hard interaction is identified using information on the directions of flight of the photons as 
determined using the longitudinal segmentation of the ECAE (calorimeter pointing), the parameters of 
the beam spot and the properties of the tracks associated with each reconstructed vertex [14]. This proce¬ 
dure ensures good diphoton mass resolution despite the presence of pileup. Pileup refers to the additional 
pp interactions that occur in the same bunch crossing as the pp ^ yy + X interaction or in the adjacent 
crossings. 

A calorimetric isolation observable is defined for photon candidates. This observable provides further 
discrimination between genuine photons and fake candidates from jets (y). It is used for two purposes 
in this analysis: to quantify the composition of the data sample in terms of true and fake photon candi¬ 
dates, and to suppress the reducible background by requiring both photons to be isolated. The isolation 
energy E™ is defined as the sum of the transverse energies of cells inside clusters of variable size (also 
called “topo-clusters” [15]) inside a cone of AE = + (Acp)^ = 0.4 around the photon candidate, 

excluding a 5 x 7 grid of EM cells in the center of the cone. 

The measurement of E^° is corrected for the following two effects: i) soft energy contributions from the 
underlying event and pileup that fall inside the isolation cone and ii) leakage of photon energy outside the 
central core, causing the isolation energy to grow as a function of photon E^. The expected soft contri¬ 
butions are subtracted from E™ on an event-by-event basis; the correction is based on a measurement of 
the properties of low-pj jets [16] in the event. The average expected contributions to Ej° from leakage 
are subtracted based on simulations. The fluctuations (from one photon shower to another) around the 
average expected leakage contribution to E^° become more important as E^ increases. In particular, the 
tail at large Ej° becomes more pronounced and would cause significant inefficiencies at large E^ if a fixed 
upper limit on Ej° were imposed in the event selection. To avoid inefficiencies, a negative E!^-dependent 
offset is added to Ej°, and the requirement E^° < 8 GeV is imposed on both photon candidates after 
the addition of the offset. In the following, the symbol E^° refers to the isolation energy with the offset 
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applied. The -dependence of the offset is chosen such that the efficiency of the requirement is 
independent of in simulated samples, and the size of the offset varies from zero at low EJ^ to -5 GeV 
SitEl^ - \ TeV. The corresponding efficiency of fhe isolafion requiremenf is 96% per phofon. For fhe pur¬ 
pose of fhe deferminafion of the background composition, a less stringent requirement of E]^° < 14 GeV 
is used. Examples of the E]^° distributions will be shown after the discussion of the backgrounds. 


4 Background estimate 


An important contribution to the background arises from prompt yy production via Standard Model pro¬ 
cesses. This contribution is irreducible and, as quantified below, represenfs the dominant source of back¬ 
ground. Another significant contribution is due to events in which one or both of the photon candidates 
arise from other objects, such as misidentified jets or electrons. This background component is domi¬ 
nated by y -I- 7 and j + j events with one or two jets reconstructed as photons. Backgrounds with electrons 
misidentified as photons - e.g. elecfron-posifron evenfs from Drell-Yan producfion, in IT/Z -i- y or in 
ft evenfs - were verified fo be negligible affer fhe evenf selecfion. 

The irreducible background is esfimafed using Monfe Carlo (MC) simulafions normalized fo fhe dafa 
in a confrol region defined below, and fhe reducible background is esfimafed using confrol samples in 
dafa. A combination of fhe leading-order (LO) evenf generator Pythia [17] and fhe nexf-fo-leading-order 
(NLO) parfon-level simulafion Diphox [18] is used fo simulafe fhe irreducible background. Monfe Carlo 
samples are generated using Pythia 8.163 wifh the CTEQ6L1 [19] parton distribution functions (PDEs) 
and the AU2 Pythia parameter tune [20]. The detector response is simulated using Geant 4 [21, 22], 
including pileup conditions similar to those observed in data. The simulated events are reweighted such 
that the number of pp interactions per bunch crossing has the same distribution as in data. The NEO cor¬ 
rection and the contributions of the fragmentation process are evaluated using Diphox 1.3.2 with the 
CTEQ6.6M [23] PDE set and are used to scale the Pythia prediction by a factor that depends on the 
generated mass of the diphoton system. The so-called box contribution gg yy through a quark loop is 
included in both the Pythia and Diphox predictions. Prom the point of view of power counting, this dia¬ 
gram is a next-to-next-to-leading-order (NNEO), i.e. contribution. But given the large gluon 

luminosity at the EHC compared to the quark-antiquark one, the contribution of the box diagram is com¬ 
parable to that of the qq yy process, which corresponds to a EO, i.e. 0{a\j^), diagram. 

To reduce the impact of the theory uncertainties in the absolute normalization of the production cross- 
section of the irreducible background, the background estimates are scaled to the data at low mass. In 
the mass region tUyy < 600 GeV, where tUyy denotes the measured invariant mass of the diphoton system, 
the most sensitive searches for narrow resonances in the diphoton channel rely on background estimates 
that make use of low-mass and high-mass sidebands [24]. The search described in the present article fo¬ 
cuses on higher diphoton masses where the sidebands, especially at high mass, provide insufficient event 
yields. This search therefore uses a different approach to obtain the background prediction at high niyy. 
Specifically, fhe simulated irreducible background is normalized fo fhe dafa in a low-mass confrol region 
affer subfracfion of fhe reducible backgrounds. This confrol region is defined fo coincide wifh fhe firsf 
22 bins in Pig. 3, i.e. as 179 < niyy < 409 GeV. To determine fhe composition of fhe dafa sample in fhe 
low-mass confrol region in ferms of irreducible and reducible backgrounds, a femplafe fif fo fhe E}^° dis- 
fribufions of fhe leading and subleading phofon candidafes is performed. A similar technique was used 
in earlier analyses of low-mass diphofon evenfs [25, 26]. The disfribufions, along wifh fhe resulf of 
fhe femplafe fif, are shown in Pig. 1 (fop and middle). The full evenf selecfion described in Sec. 3, wifh 
fhe relaxed requiremenf E™ < 14 GeV, has been applied. The fif is performed simulfaneously fo fhe 
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two distributions, and the normalizations for four background components are allowed to float in the fit, 
namely the irreducible yy component plus three reducible components: y + j (subleading photon candi¬ 
date is fake), j + y and j + j. Templates for the distribution of true photons and of fake photons are 
determined from data. The template for fake photons is obtained using a sample of photon candidates 
that pass the non-tight selection: photon candidates must satisfy the loose identification criteria and fail 
to meet a subset of the criteria in the tight selection. This subset has been chosen to include criteria that 
have been found to be only weakly correlated with Ej°. The template for true photons is obtained from 
the sample of tight photon candidates after subtraction of the fake-photon template normalized to the tight 
data sample at Ei^° >10 GeV where fake photons dominate. 

The fit is performed separately in five subsamples of the data, because the relative contributions of the 
four background components vary significantly from one subsample to another. The five subsamples cor¬ 
respond to different combinations of pseudorapidities of the leading and subleading photon candidates: 
both photons in the “central” region (I 77 I < 1.37), leading photon central and subleading photon in the 
“endcap” region (1.52 < It/I < 2.37), leading photon in endcap and subleading central, both photons in 
same endcap and both photons in endcaps of opposite-sign pseudorapidity. In Fig. 1, the data from the 
five subsamples, as well as the corresponding fit results are combined, and only the sum of the photon 
and jet components over the four background components is shown. To illustrate the broadening of the 
F'™ distribution at large E^ described above. Fig. 1 also includes the Ei^° distribution for photons with 
E\ > 400 GeV. 

The backgrounds in the search region, i.e. at masses higher than the low-mass control region, are pre¬ 
dicted using estimates of the tUyy shape of each background component and the normalizations obtained 
in the low-mass control region described above. As discussed above, the shape of the irreducible compo¬ 
nent is obtained using simulation. The shapes of the reducible components are obtained from data control 
samples. The control sample used to extract the shape of the y + j background is selected in the same 
way as the signal sample, except that the subleading photon candidate has to satisfy the loose identifica¬ 
tion criteria and fail to meet the tight criteria. The control samples for the j + y and j + j backgrounds 
are defined using the same approach. Since these data control samples contain relatively few events at 
high niyy, a fit to a smooth function of the form f{myy) = p\X where the pi are free 

parameters, is used to extrapolate the shapes of the reducible backgrounds to high niyy. This functional 
form describes well the shapes of the data in the control samples of this analysis. It was also successful 
in describing the shapes of the background control samples in earlier searches in the dijet [27, 28], dilep¬ 
ton [29], photon-jet [30] and diphoton [5] channels. 

The shapes of the distribution of the isolation energy for genuine photons extracted from data (in bins of 
Ely) are compared to the shapes predicted using simulation. The impact of the small observed differences 
in the isolation distributions on the predicted niyy shape of the irreducible background and on the signal 
efficiency is quantified and taken into account in the systematic uncertainties. 


The uncertainties in the prediction of the total background are shown in Fig. 2 as a function of niyy. The 
uncertainties in the shape of the niyy distribution of the irreducible background are dominated by the 
uncertainties in the PDF set used for the Diphox simulation. The set of PDF variations provided with 
the CTEQ 6 . 6 M PDF set are used to propagate the PDF uncertainty to the shape of the irreducible back¬ 
ground. The comparatively small difference between the irreducible background shapes obtained with the 
CTEQ 6 . 6 M and the MSTW2008 [31] PDE sets is added in quadrature to the uncertainty. An additional 
uncertainty arises from higher-order contributions that are not included in the Diphox generator. This 
uncertainty is evaluated by varying in a correlated and anticorrelated way the renormalization scale and 
the factorization scales in the hard scattering process and in the photon fragmentation function [18] by a 
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Figure 1; (Top) Distribution of the isolation energy is™ for the leading photon candidate in events in the low-mass 
control region. (Middle) Corresponding distribution for the subleading photon. The black points represent the 
data. The blue line shows the result of the template fit to the data. The red and green lines show the templates for 
photon and jet background respectively, normalized according to the fit result. The red and green bands represent 
the uncertainties in the photon and jet contributions due to the uncertainties in the template shapes. (Bottom) 
Distribution of F™ for photon candidates with E\ > 400 GeV at any diphoton mass. The F™ values always 
include the F^-dependent offset discussed in the text. 
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Figure 2: Uncertainty on the total background estimate as a function of niyy. The individual contributions due to 
the uncertainties in the irreducible and reducible components, as well as the contribution from the uncertainties in 
the normalization of the different background components (labeled “background composition”), are also shown. 


factor two around their nominal value, which is set to the mass of the diphoton system. A small uncer¬ 
tainty arises from residual imperfections in the simulation of the isolation energy E^^°. The uncertainties 
in the shape of the reducible background arise from the finite size of the corresponding control samples, 
and from the extrapolation of the background shapes from the control samples to the signal region. The 
latter are assessed by varying the definition of the loose selection requirement that is used to define the 
control samples. The uncertainties in the normalizations of the individual background components are 
dominated by the uncertainties in the template shapes that are used in the fit to the isolation energy. They 
are assessed by varying the definition of the non-tight selection requirements that are used to obtain the 
template shapes. 


5 Signal response 


Simulated event samples of RS graviton signals at different values of mo* are used to study the response of 
the detector for signal events. The samples are generated using Pythia 8.163 with the CTEQ6L1 PDF set, 
the AU2 tune and the same detector and pileup simulation as for the SM yy samples described above. The 
sources of uncertainty in the expected signal yield are summarized in Table 2. The systematic uncertainty 
is dominated by the uncertainty in the efficiency due to residual imperfections in the simulation of the 
observables that are used in the tight photon identification criteria and in the simulation of the isolation 
energy. Smaller contributions arise due to the limited size of the simulated samples and the determination 
of the trigger efficiency. The uncertainty in the integrated luminosity is also propagated to the expected 
yield. Uncertainties in the RS resonance shape due to our current knowledge of the photon energy scale 
and resolution as well as imperfections in the simulation of the pileup conditions were verified to have a 
negligible effect on the result. 






Table 2; Summary of systematic uncertainties on the expected signal yield, excluding those on the production cross 
section and the acceptance. The total systematic uncertainty is obtained by summing the individual contributions 
in quadrature. 


Source 

Uncertainty in signal yield [%] 

Integrated luminosity 

2.8 

MC statistics 

1.0 

Trigger efficiency 

1.0 

Photon ID efficiency 

3.0 

Photon isolation efficiency 

0.3-2.1 

(for mo* - 500-3000 GeV) 

Total 

w5 



Figure 3; Observed invariant mass distribution of the selected diphoton events (black dots; the vertical and hori¬ 
zontal axes on logarithmic scales). To compensate for the rapid decrease of the spectrum, the bins are chosen to 
have constant logarithmic width. Specifically, the ratio of the upper to lower bin boundary is equal to 1.038 for 
all bins, and the hrst bin starts at 179 GeV. Superimposed are the SM background prediction including irreducible 
and reducible components and two examples of signal predictions. The low-mass control region is indicated by an 
atTow. (Bottom) Bin-by-bin significance of the difference between data and background expectation. 
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Table 3: Number of events expected from the reducible and irreducible background components as well as the total 
background and observed number of events in different mass bins. Each bin boundary in the table corresponds to a 
bin boundary in Fig. 3, but the binning of the table is more coarse than in the figure. Since the combined statistical 
and systematic uncertainties are strongly anticorrelated between the two components, the uncertainty in the total 
background prediction is smaller than the sum in quadrature of the uncertainties of the individual components. 


Mass window 

Background expectation (number of events) 

Observed 

[GeV] 

Irreducible 

Reducible 

Total 

events 

[179,409] 
Control region 

23800 ± 2400 

9100 ±2400 

32866 

32866 

[409,513] 

1070 ± 110 

400 ± 100 

1463 ± 27 

1465 

[513,596] 

369 ± 37 

129 ± 34 

498 ± 12 

524 

[596,719] 

240 ± 24 

74 ±20 

314.4 ± 8.8 

335 

[719,805] 

75.8 ±7.7 

20.6 ± 5.5 

96.4 ± 3.2 

99 

[805,901] 

46.6 ± 4.8 

11.5 ±3.2 

58.1 ± 2.1 

60 

[901,1009] 

28.2 ± 3.0 

6.3 ± 1.8 

34.5 ± 1.5 

33 

[1009,1129] 

16.8 ± 1.9 

3.4 ± 1.0 

20.2 ± 1.0 

15 

[1129,1217] 

6.92 ± 0.89 

1.35 ± 0.46 

8.27 ± 0.45 

7 

[1217,1312] 

4.85 ± 0.73 

0.88 ± 0.39 

5.74 ± 0.36 

3 

[1312,1415] 

3.11 ±0.54 

0.58 ± 0.28 

3.69 ± 0.25 

0 

[1415,1644] 

3.39 ± 0.59 

0.61 ± 0.29 

4.00 ± 0.32 

5 

[1644,3000] 

2.12 ±0.61 

0.41 ± 0.22 

2.52 ± 0.30 

3 


6 Results 

The observed diphoton mass spectrum is shown in Fig. 3, together with the background expectation and 
the predicted signal for two examples of values of the RS model parameters. The acceptance (i.e. the frac¬ 
tion of simulated events in which both photons satisfy the requirements on It/I and described above) 
for the RS scenario with mo* - 2.0 TeV and kjMpx - 0.1 is 89.3% and the selection efficiency for 
events within acceptance is (55.9 ± 2.6)%. The bottom inset in Fig. 3 shows, bin-by-bin, the statistical 
significance of the difference between the data and the expected background. This significance is cal¬ 
culated using the prescription in Ref. [32] and is plotted as positive (negative) for bins with an excess 
(deficit) of data. The observed event yields and expected backgrounds are also summarized in Table 3. 
No significant deviation of the data from the expected background is observed. To further quantify the 
level of agreement, the BumpHunter [33] algorithm is used to search for any “bump” on top of the ex¬ 
pected background, i.e. any excess of events localized in a relatively narrow mass region. Specifically, 
the BumpHunter algorithm scans the mass distribution in windows of progressively increasing width in 
order to identify the window with the most significant excess of observed events over the background 
anywhere in a predefined search region. In this analysis, the binned myy spectrum from Fig. 3 is used, 
the search region is defined as 409 < myy < 3000 GeV, and the window size is allowed to vary from one 
bin to half the number of bins in the search region. The most significant excess is found in the region 
494 < myy < 745 GeV. The probability to observe anywhere in the search region an excess that is at least 
as significant as the one observed in data and that arises from fluctuations in the background is 58%. This 
confirms the absence of any significant excess in the data. 
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Table 4: Expected and observed lower limits at 95% CL on me* for different values of klMp\. 


kj Mp\ 

-2cr 

Expec 

-Icr 

:ted limit 

central 

[TeV] 

-l-lcr 

-l-2cr 

Observed 
limit [TeV] 

0.010 

1.45 

1.38 

1.30 

1.20 

1.11 

1.41 

0.020 

1.80 

1.78 

1.70 

1.60 

1.49 

1.62 

0.030 

2.02 

2.01 

1.96 

1.86 

1.77 

1.92 

0.040 

2.17 

2.16 

2.14 

2.04 

1.94 

2.15 

0.050 

2.29 

2.28 

2.27 

2.17 

2.08 

2.28 

0.060 

2.38 

2.38 

2.37 

2.27 

2.19 

2.38 

0.070 

2.47 

2.46 

2.45 

2.37 

2.27 

2.46 

0.080 

2.54 

2.54 

2.53 

2.46 

2.38 

2.54 

0.090 

2.61 

2.60 

2.60 

2.53 

2.45 

2.60 

0.100 

2.67 

2.66 

2.66 

2.60 

2.51 

2.66 


As no evidence for a signal is found, upper limits at 95% confidence level (CL) are set on the production 
cross section times branching ratio, crxBR(G* —> yy), for the lightest RS graviton. A Bayesian approach 
is used to calculate the limits, and the corresponding computer program is implemented using the toolkit 
from Refs. [34, 35] . A likelihood function is used to describe the probability to obtain the rtiyy distribution 
observed in data, given the background and signal predictions. The likelihood is defined as the product 
of the per-bin probabilities over all bins in the search region {myy > 409 GeV) in Fig. 3. In each bin, the 
probability is computed from a Poisson distribution given the observed number of events and the model 
prediction, where the model prediction is the sum of the expected numbers of background and signal 
events. The expected signal yield depends on me* and cr x BR{G* —> yy). Using Bayes’ theorem, the 
likelihood is converted into a posterior probability. A uniform positive prior in cr x BR{G* —> yy) is 
assumed, and systematic uncertainties are incorporated using nuisance parameters with Gaussian priors. 
For a given value of , the posterior probability is reduced to one parameter of interest, cr x BR(G* 
yy), and the 95% CL limit on cr x BR{G* —> yy) is obtained via integration of the posterior probability. 
The resulting limits on cr x BR{G* —> yy), as a function of the mass hypothesis (value of mo*), are shown 
in Fig. 4. The phenomenology of the RS model can be described in terms of the parameters me* and 
kjMpi, i.e. cr X BR(G* yy) can be calculated [36, 37] for a given pair of values for these parameters. 
The results of this calculation are represented by the lines for given values of k/Mpi in Fig. 4. They are 
obtained using the Pythia generator, which implements the calculations from Ref. [37]. A A'-factor is 
applied to these predictions to account for QCD corrections at NLO. The value of the A'-factor varies 
between 1.6 and 1.9, depending on mo* and k/Mpi. These results for the A'-factor were provided by 
the authors of Refs. [38, 39], using updated calculations at = 8 TeV. Given these theory predictions 
for cr X BR{G* —> yy), the upper limits on cr x BR{G* —> yy) from Fig. 4 can be translated into upper 
limits on kjMpx. The results are shown, as a function of tuq-, in Fig. 5. Alternatively, the upper limits 
on cr X BR{G* yy) can be translated into lower limits on the graviton mass. The results are shown 
in Table 4, for different values of k/Mpp The limits presented in Fig. 5 and in Table 4 do not take into 
account any theory uncertainties in the prediction of cr x BR{G* —> yy) for the RS signal. The size of 
one of the dominant theory uncertainties (the PDF uncertainty) is included in Fig. 4; uncertainties due to 
missing higher-order contributions in the calculation of the A'-factor are estimated to be of comparable 
size. 
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[TeV] 


Figure 4: Expected and observed upper limits on cr x BR{G* —> jy) expressed at 95% CL, as a function of the 
graviton mass. At large me-, the -Icr and -2cr variations of the expected limit tend to be particularly close to the 
expected limit. This is expected, since signals with high rriQ- would appear in regions of niyy where the background 
expectation is small and the Poissonian fluctuations around the mean expected background are highly asymmetric. 
The curves show the RS model prediction for given values of k/Mpi as a function of mg.. They are obtained using 
the Pythia generator plus a /T-factor to account for NLO corrections (see text). The thickness of the theory curve 
for k/Mpi = 0.1 illustrates the PDF uncertainties expressed at 90% CL. 


7 Conclusions 

In summary, a search for high-mass diphoton resonances has been performed using the full 2012 dataset 
collected by the ATLAS detector at the LHC (20.3 fb“^ of pp collisions at = 8 TeV). No significant 
excess over the expected background is observed, and upper limits on the production cross section times 
branching fraction cr x BR{G* —> yy) for narrow resonances are reported as a function of the resonance 
mass. High-mass diphoton resonances are predicted e.g. in models that postulate the existence of extra 
spatial dimensions in order to address the hierarchy problem, and could be observed at the LHC. Limits on 
the mass of the lightest graviton state in the framework of the Randall-Sundrum model of extra dimensions 
are derived from the limits on cr x BR{G* —> yy). The only free parameters of this model are the mass 
of the lightest graviton {mQ>) and the coupling to the SM fields (klMpi), i.e. cr x BR{G* yy) can be 
calculated once fhe values of fhese fwo paramefers are specified. The limifs can fherefore be expressed 
in terms of mo* for a given value of kIMpi. A lower limif of 2.66 (1.41) TeV af 95% confidence level is 
obfained on fhe mass of fhe lighfesf graviton for a coupling kjMpx - 0.1 (0.01). The resulfs reported in 
fhis article consfifufe a significanl improvemenf over fhe resulfs in Ref. [5] obfained at ^fs -1 TeV. The 
upper limifs on cr x BR{G* —> yy) are reduced by a faefor of 2.1 (3.6) for a resonance wifh a mass of 
0.5 TeV (2.75 TeV). Furfhermore, fhe new limifs are expressed in ferms of fhe production cross section 
af = 8 TeV, compared to = 7 TeV for fhe old limifs, and fhe production cross section for any 
new heavy particle is expected fo be larger af = 8 TeV. The lower limifs on fhe mass of fhe lighfesf 
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Figure 5; Expected and observed upper limits on kjMpi expressed at 95% CL, as a function of the graviton mass. 


RS graviton are also significantly improved compared to the previous analysis. 
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